Quick and dirty methods for studying black-hole resonances. 
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We discuss simple integration methods for the calculation of rotating black hole scattering res- 
onances both in the complex frequency plane (quasinormal modes) and the complex angular mo- 
mentum plane (Regge poles). Our numerical schemes are based on variations of "phase-amplitude" 
methods. In particular, we discuss the Priifer transformation, where the original (frequency domain) 
Teukolsky wave equation is replaced by a pair of first-order non-linear equations governing the intro- 
duced phase functions. Numerical integration of these equations, performed along the real r,-axis 
(where r, denotes the usual tortoise radial coordinate), or along rotated contours in the complex r,- 
plane, provides the required iS-matrix element (the ratio of amplitudes of the outgoing and ingoing 
waves at infinity). Miiller's algorithm is then employed to conduct searches in the complex plane 
for the poles of this quantity (which are, by definition, the desired resonances). We have tested 
this method by verifying known results for Schwarzschild quasinormal modes and Regge poles, and 
provide new results for the Kerr black hole problem. We also describe a new method for estimating 
^ the "excitation coefficients" for quasinormal modes. The method is applied to scalar waves moving 

'^'n , in the Kerr geometry, and the obtained results shed light on the long-lived quasinormal modes that 

■ exist for black holes rotating near the extreme Kerr limit. 
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During the last few decades it has become clear that resonant phenomena play an important role in black hole 
physics Numerical experiments, pioneered by Vishveshwara more than thirty years ago, and continued since 
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then by many other authors, clearly show the dominance of the so-called quasinormal mode (QNM) oscillations in 
, most dynamical processes involving black holes. The importance of this phenomenon is further emphasised by the 
■ possibility that future observations of QNM "ringing" following, say, the merger of two black holes by the network of 
' interferometric gravitational wave detectors may provide direct identification of black holes. This exciting prospect 

has motivated a large number of investigations into the nature of the QNMs. An exhaustive list of papers on the 
^ ' subject can be found in Refs. 

(^JT), The spectrum of QNM frequencies has been studied extensively in the literature, especially for the Schwarzschild 
black hole. It is well-known that the actual mode calculation is a non-trivial task. Recall that any physically acceptable 
, solution to the relevant (homogeneous) wave equation takes the form of a purely "ingoing" wave ~ exp(— zcjr*) 
• (where r» is the tortoise radial coordinate) at the event horizon and a mixture of "ingoing/outgoing" waves 
|h [ exp(— iwr*) -t- Aout exp(ia;r^,) far away from the black hole (such solutions exist because the effective curvature 

potential falls off asymptotically towards both the horizon and infinity, i.e. when r, — > ±oo). A QNM solution 
corresponds to purely outgoing waves at infinity. In other words, to search for a QNM one must first calculate the 
asymptotic amplitudes ^out and ^in, and then try to find a frequency that satisfies the required QNM boundary 
condition = 0. However, direct numerical integration of the wave equation leads to serious numerical problems. 
Since the QNMs are expected to be damped due to gravitational-wave emission the corresponding solutions must 
grow exponentially as one approaches infinity and the horizon (on a i-constant hypersurface). In order to identify a 
QNM one must be able to filter out the solution that decays exponentially at infinity. In other words, the numerical 
calculation requires exponential precision. 
_ Since the first "direct" calculation of the slowest damped Schwarszchild QNMs by Chandrasekhar and Detweiler 
various alternative schemes have been devised. Detweiler |^ was able to determine the fundamental Kerr QNM 
with adequate accuracy by examining the functional form of the asymptotic amplitudes along the real frequency axis. 
Ferrari and Mashhoon *?! exploited the correspondence between resonant and bound states of the inverted potential 
to obtain approximate black hole modes. WKB techniques were successfully used to calculate the first few least 
damped modes by Schutz and Will Q and several subsequent authors Jj^ 0, 0, ^ . High precision QNMs for both 
Schwarzschild and Kerr black holes were first calculated by Leaver |l3lll4| | by means of a continued fraction approach. 
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* Presently performing national service with the Greek army. 
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His results for the Schwarzschild case were later verified (and extended) by Andersson [T^ via the phase-amplitude 
formalism and NoUert and Schmidt (161] who used a Laplace transform/analytic continuation method. Moreover, the 
semi-analytic phase-integral formalism jl^ has proved to be a valuable tool for the calculation and interpretation of 
QNMs [la. 

More recently, Onozawa extended Leaver's calculation of the QNMs of rotating black holes and studied the 
behaviour of the mode frequencies as the rotation parameter is varied. The "long-lived" QNMs of a near extreme 
Kerr black hole were discussed by the present authors '2l'| , using analytic expressions for the asymptotic amplitudes 
derived by Press and Teukolsky 22]. That study complemented a calculation by Detweiler who was the first to use 
these amplitudes to derive a QNM-condition for near extreme rotation, and moreover, solved it analytically for the 
extreme case . Finally, Berti and Kokkotas ^3 have used the continued fraction approach to study the behaviour 
of the high overtone QNMs of a Kerr black hole. Their study unveils a complicated, and quite puzzling, structure. 

In the last few months there has been renewed interest in the spectrum of black hole QNMs. Motivated by the 
suggestion that the asymptotic behaviour of the high-overtone QNM freauenci es .24..25J may be related to the so-called 
Barbero-Immirzi parameter which is relevant for loop quantum gravity |2fil l27|]. there has been several attempts to 
approximate the high-frequency part of the spectrum. Some progress has been achieved by linking the high-overtone 
spectrum with the black hole's quantum properties [28i] . Although we will not discuss this issue further in the present 
article, it is an interesting problem that need to be investigated in complete detail. 

The notion of Regge poles has (so far) hardly at all been discussed in the context of black holes. The Rcgge poles 
are close relatives to the QNMs 29] , and correspond to those complex values of the angular momentum / for which 
the wavefunction is purely outgoing at infinity and purely ingoing at the horizon. The Regge pole calculation poses 
an easier problem than that of finding QNMs since the corresponding wavefunction is well behaved at infinity (and 
the horizon) for complex / (as long as the frequency is kept real) . In quantum theory, the complex angular momentum 
(CAM) formalism (in which Regge poles play a central role) provides an alternative (and usually more transparent) 
description of scattering problems (Sfl]. Some years ago, one of us apphed the CAM theory to the study of scattering 
of massless scalar waves by a Schwarzschild black hole . In that study, Regge poles and the corresponding 

residues were calculated via the phase-integral method. The CAM theory has also been used by Chandrasekhar and 
Ferrari |33j] in an investigation of axisymmetric pulsations of a relativistic star. 

In this work we discuss numerical methods for calculating both QNMs and Regge poles. The key motivation for 
our study was the desire to devise methods that were easy to implement (hence "quick and dirty") and yet able to 
produce numerically accurate results. To achieve this goal we express the Teukolsky wave function in terms of suitable 
phase- functions (see Section III below). The integration of the resulting first-order equations readily yields the desired 
ratio of asymptotic amplitudes. In order to retain sufficient numerical precision to be able to identify the ingoing 
wave at infinity for complex frequencies, while ensuring that the phase-functions remain smooth and well behaved, the 
integration is generally performed along straight lines in the complex r^, plane. The slope of the integration contour 
is chosen in such a way that the relative magnitude of the out- and ingoing wave solutions is preserved both near 
the horizon and at infinity. The end result is a multi-purpose scheme that (with minor modifications) facilitates the 
calculation of QNMs and Regge poles as well as scattering phase-shifts 's^. 

The remainder of the Paper is organized as follows: In Section II we briefly discuss the formal deflnition and 
physical interpretation of black hole resonances. Our numerical approaches to the problem are described in full detail 
in Section III. Section IV contains numerical results for QNMs (Section IVA) and Regge poles (Section IVB). Section V 
provides a summary and some final comments. Geometric units G — c — 1 are adopted throughout the paper. 



II. BLACK HOLE RESONANCES: DEFINITION AND PHYSICAL MEANING 

The notion of scattering resonances is probably most familiar from quantum theory |35|. The fact that it is also 
relevant in black hole physics (in terms of QNMs and Regge poles) is not particularly surprising given that linear 
perturbations of black hole spacetimes are governed by Schrodinger type equations. The corresponding effective 
potential barriers are a combination of the black hole's gravitational pull and a repulsive centrifugal term ,36]. This 
means that one can formulate a scattering problem for black holes and investigate phenomena analogous to those in 
the quantum case. For simplicity we will mainly consider the case of scalar waves in this paper. This is not a severe 
restriction, however, since more "realistic" cases, eg. gravitational perturbations, can be dealt with in exactly the 
same way. 

After adopting the Boyer-Lindquist coordinate frame, and separating the angular dependence in the standard way, 
we arrive at the Teukolsky equation )37j. which governs the radial component 4'{r) of a scalar perturbation in a Kerr 
black hole spacetime, 

-7-Y + Qirn[r*,uj,ajip = (1) 
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We have also assumed a harmonic time dependence e The tortoise coordinate is defined as 

dr, _ + 
dr " A 

where A = (r — t-i. ) (r — r_ ) and r± = M ± {M^ — a^)^/^ denote the black hole's inner (Cauchy) and outer (event) 
horizon. Explicitly we have 

^■In — -1 ^In — -1 (3) 



The effective potential Qim is given by 

^ , , + i2amLu - a^uj^ - Ei„^)A dG 

Qi^n{r.,-,-) = J^^T^2 rf^T-G (4) 

with X = (r^ + a^)uj — am and G = rA/ (r^ + a^)^. It has the following asymptotic behaviour 

J /c^ as r -> r+ (r, ^ -oo) , . . 

\ w2 as r ^ +C30 (r, oo) , 

where k — to — muj^ with = a/2Mr+ denotin g th e angular frequency of the event horizon. The eigenvalue of the 
relevant angular equation is denoted by Ei^m cf. j33| . 

As a result of |(SJ|, the solution of will (to leading order) be a linear combination of simple exponentials as 
r* ±oo. A physically acceptable solution describes a wave propagating towards the horizon at — oo, i.e. a 

purely "ingoing" solution. We can normalize this solution according to 

^ \ Aine-^'^'-- + ^oute*"'^* as ^ +00 . 
Another useful solution is purely "outgoing" at infinity, 

The Wronskian of any two linearly independent solutions of is constant. Using the two solutions above we easily 
get 

l^(^»,V"^) = 2ic^A„ (8) 
Finally, we define the tS-matrix, a quantity familiar from scattering theory [35j |. as the ratio 

/i out 

S{u;,l)^i-lY+'^=cM2^Sl) (9) 

where Si defines the scattering phase-shifts. 

One can define two classes of "resonances" that correspond to purely ingoing waves at the horizon and purely 
outgoing waves at infinity. Loosely speaking, such solutions are "characteristic" of the black hole spacetime since they 
do not depend on waves coming in from infinity. The first class of resonances are the QNMs. They correspond to 
complex frequency poles ujn of S (for integer I). The QNM frequencies (obviously) correspond to solutions such that 
the Wronskian ((SJ vanishes. The second possible class of resonances are the Regge poles. To identify these one allows 
the angular momentum (I) to assume complex values (while keeping lo real) and searches for poles /„ of S. The close 
connection between these classes of solutions has recently been discussed by Decanini, Folacci and Jensen [29| . 

It is well-known that the numerical determination of QNMs requires considerable care. Since we are discussing 
(stable) damped modes, the QNM frequencies will have negative imaginary parts. It follows then, that each QNM 
eigenfunction e~*""*V''°('^n, ?"*) will grow exponentially both towards infinity and at the horizon. Since one must be 
able to filter out the ingoing-wave solution at infinity in order to identify a QNM, the calculation requires exponential 
precision. It should be pointed out that this problem is not physical: It arises because we are viewing only a part of 
the overall problem, as is clearly shown in discussions of the initial- value problem for perturbed black holes |3llla. l40| . 
Those studies reveal the true physical meaning of QNMs: They describe part of the late time behaviour of the 
perturbing field as observed at some fixed distance from the black hole They can be viewed as arising as a part 
of the field is temporarily "trapped" near the scattering center (the black hole) and leaks out to infinity. The same 
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picture arises in the study of QNMs of much simpler systems such as static objects scattering waves in flat spacetime 

The QNMs arise in scattering situations that can be thought of as "time dependent" (strictly speaking, the evolution 
and scattering of initial data with compact support). Alternatively, one might consider a time- independent scattering 
"experiment" . Namely, the scattering of monochromatic plane waves (where the time dependence is trivially given by 
~ exp(— zciji) for any given lu). In this case, the most important "observable" quantity is the differential cross section, 
which is a measure of the scatterer's size as seen by waves scattered into a certain solid angle |35j| . The field at infinity 
consists of the incoming plane wave plus a scattered piece, and can be written 

V'(n -^+00)^V'planc + /(e)^ (10) 

where we have assumed, for simplicity, axisymmetric scattering (for a treatment of non-axisymmetric scattering 
by rotating black holes, see [s^). The cross section can be expressed in terms of the (complex- valued) scattering 
amplitude f{9), 

s ^ mi' (11) 

For scalar waves scattered off a Schwarzschild black hole one finds that |41|| , 

^ oo 

fie) = 77-Y.^^l + '^){Si-micos9) (12) 

from which we see that all the relevant scattering information is encoded in the 5-matrix. 

The angular scattering problem is usually studied in terms of the phase-shifts Si. However, it is not unusual for 
the sum in (|12l) to be slowly converging (black-hole scattering being a example, although this difficulty can to some 
extent be avoided following the prescription in ^3)- many situations it is therefore desirable to have a different 
representation for the scattering amplitude. The most celebrated alternative is the CAM approach "sd]. By means 
of a Watson-Sommerfeld transformation the sum over / in (|12|l can be rewritten as an integral along a contour in the 
complex /-plane This integral naturally splits into two terms, 

f{0)^ fp{e) + fBio) (13) 

The first term is a sum over the residues of the poles In of S; the Regge poles [s^- It is explicitly given by, 

/pW-^-E^^^/'a-os.) (14) 
u! ^-^ cos7r(/„ + 1/2) 

n— 

where r„ is the residue of each Regge pole, i.e. is defined by 

(15) 

in the vicinity of the n**^ pole, and P; is a Legendre function of complex degree. The second term in (|13|) is a "back- 
ground integral" , which is known to be less relevant as far as diffraction effects are concerned. As was demonstrated 
in [s^l, only the first few poles are required for determining the scattering cross section for angles away from the 
forward direction (the direction of incoming plane wave). Despite the fact that the problem of scattering by rotating 
black holes has not yet been studied within the CAM framework, it is natural to define Kerr black hole Regge poles 
as singularities of S in the complex I plane while keeping uj real and m an integer. 

The physical interpretation of Regge states is discussed in standard textbooks [s^, 113, cf. also the discussion 
of scattering of plane waves by a Schwarzschild black hole in |32j |. By approximating the Legendre function for 
|/„ -|- 1/2| ^ 1 and 6 not too close to or tt one can show that 

-,1/2 
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2tt{1„ + l/2)sin6' 



|gj(i„ + l/2)(7r-e)-j7r/4 _^ g-j(i„ + l/2)(7r-e)+j7r/4| j^^g^ 



From this formula it follows that each Regge pole can be interpreted as a combination of two "surface" waves traveling 
around the black hole in opposite directions. These waves decay exponentially as they propagate, with an "angular 
life" ~ l/Im(/„ + 1/2). The real part of the Regge pole is associated with the distance i?„ from the black hole at 
which the angular decay occurs, Re (/„ -|- 1/2) « cuRn (note that this "localization principle" is formally valid only for 
I 3> 1 113 )• As can be verified using the numerical results presented in Section VI, the value of i?„ lies very close to 
the critical impact parameter for null geodesies (which is be = ^VSM for a Schwarzschild black hole). This suggests 
that (as one would have expected) Regge states describe waves temporarily trapped in the vicinity of the unstable 
photon orbit, leaking to infinity and to the horizon, as they travel around the black hole. 
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III. TWO "QUICK AND DIRTY" SCHEMES 



In this Section we will discuss two closely related methods for determining black-hole resonances. Both methods 
are based on the idea of integrating slowly varying phase-functions rather than the, possibly rapidly, oscillating 
wavefunction itself. The representation of the solutions to Schrodinger type equations in terms of phase-functions was 
first discussed by Milne and, a long time later, by Newman and Thorson 01- The second key ingredient in our 
schemes is the use of complex values of the radial coordinate. By analytically continuing the resonant solutions into 
the complex coordinate plane one can suppress the divergences associated with (in particular) the QNM boundary 
conditions. In essence, our methods are close relatives to the phase-amplitude method that was developed by one of 
us to study of scattering resonances. The main difference is that the methods we describe here make use of much 
simpler integration contours, and are consequently much easier to implement. 



A. The Priifer method 



The family of the so-called Priifer transformations has been widely used in the numerical analysis of Sturm-Liouville 
problems 45] . It has also been applied to the calculation of resonances and phase-shifts in standard quantum scattering 
problems, see for example Refs. (46l l47j . 

Consider a Schrodinger-type equation, where either the frequency lo or the angular momentum I are allowed to take 
complex values, 

+ Q(x,Z,c^)V' = O (17) 

We consider cases where the radial coordinate x spans the entire real axis (in the black hole case x corresponds to 
the tortoise coordinate r,), and the effective potential Q is a single barrier with asymptotic behaviour 

Q ^{ 1^ as a; --oo 

^ ]^uj'^ as X — > +00 ^ ' 

where q = q{uj). It follows then, that the function ip^Xjlju) will be a linear combination of exponentials as x ^ ±oo. 
Following the discussion in the previous Section, let us assume that we are interested in a solution to (|17|l which is 
purely "ingoing" as a; — > — oo and mixed ingoing/outgoing as x — > -t-oo. That is, we want to determine the "in" 
solution, see eqn. @, which can be written 

^ r e-^^ as X -oo , 
^ I ±>sm[ujx + Q as X — > -l-oo , ^ ' 

where ( and B are some complex constants. 

We can write the exact solution of H17|l in the form 

■0(x, /, uj) = exp[ y P(x, I, uj)dx] (20) 

where P is the logarithmic derivative of -ip (a prime denotes derivative with respect to x) , 

^ = P (21) 

This is our first phase-function. Obviously, it should obey the boundary condition P — > —iq for x —^ —oo. Alterna- 
tively, we can express the function -0 and its derivative via the transformation, 

ip{xj,uj) = Bsm[ojx + P{x,l,uj)] (22) 
ip'{x,l,LLj) = Bujcos[u!x + P{x,l,oj)] (23) 

where P is the Priifer phase function (which should have C as its limiting value as x +oo). 
The two phase functions satisfy non- linear first order differential equations; 

HP 

^ + P^+Q = (24) 
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dx \ uj 
and they are directly related by 

1 \iP~uj 
P = -UJX + — log — - — 
2i iF + Lo 



'^^ + (cc>-^)sin2(a;x + P)=0 (25) 



(26) 



P = wcot(cja; + P) (27) 

A major advantage of using phase-functions can be seen immediately from Eqns. H24|) and H25|) . assuming for the 
moment that we are dealing with real lo and I. Using (|18|l . we see that (|24ll and H25|) lead to 



dP_ 
dx 



as X ^ — oo (28) 



dP 

— w as X ^ +00 (29) 
dx 

Hence, the functions P and P are slowly varying for large negative and positive values of x, respectively. This 
means that the oscillatory behaviour of the original wavefunction in those regimes [which would cause the build up 
of unnecessary numerical errors in a direct integration of (|17|l ] has been effectively removed. 
Finally, knowledge of the P phase-function at 2: — > -|-c» provides us with the iS-matrix, 

5(Z,w) = (-l)'+ie2< (30) 

Obviously, the described method can be applied to the black- hole problem with x = r.^,, q = k and (|17|l becoming 
the Teukolsky radial equation 

At the practical level, the numerical integration for the P function shoud be initiated at a point (say), close 
to the event horizon. Assuming that Q is a slowly varying function for large |r*|, which is certainly the case in the 
applications under consideration, we can initiate the integration of (|24|) from the lowest order WKB approximation 

Pw±iQi/2 as|n|^cx) (31) 

(alternative approximations, eg. a power series expansion, will work equally well). In this way we can readily find 
a solution that corresponds to the (analytically continued) ingoing wave boundary condition at the horizon. This 
solution is then integrated towards the matching point r™**. As the integration reaches the vicinity of the potential 
peak (the local minimum of Q in the Kerr case) we switch to the P function, and then carry on the integration to 
"infinity" . It is mandatory to change phase- function if one wants to ensure that the integrated function remains smooth 
and well behaved. This is due to the fact that the phase-functions undergo a drastic change in their behaviour (they 
become rapidly oscillating and may even have singularities) as they are extended through the region of the potential 
peak. This is a manifestation of the so-called Stokes phenomenon — the "switching on" of small exponentials in the 
solution to an equation of form lfT7|l . well known from WKB-and phase-integral theory (see also the discussion 

in Appendix 1 of ^3)- This situation is demonstrated in Fig. ^ where we plot the functions P{r*) and P{r*) for 
typical parameters I = m = 2, ojM = 0.5, a = 0.5M. 

The Priifer approach has already been applied in the context of scattering of plane scalar waves by rotating black 
holes 01 ■ In that study, where both lo and I were assumed to be real, the use of the phase- functions described above 
facilitated the computation of accurate high frequency phase-shifts (which would be problematic if the original wave 
equation were to be integrated |42j). We can follow exactly the same procedure when I is allowed to acquire complex 
values (as long as uj remains real). In this way we have calculated the first few Regge poles of both Schwarzschild and 
Kerr black holes. The numerical results will be discussed in the following Section. 

Complications arise when to becomes complex while I is kept real (integer). From Eqn. (|25|l we see that the term 
which represents an ingoing (outgoing) wave becomes exponentially small (large) at large x (for Im a; < 0). A similar 
situation appears close to the horizon. Therefore, we are still facing one of the problems associated with a direct 
integration of the wave equation. An efficient and elegant way to avoid this difficulty is to consider the problem in 
the complex r*-plane. The success of this approach has been repeatedly demonstrated with the application of the 
phase-integral/amplitude formalisms to the calculation of black- hole QNMs p^ll9l l55| . The main idea is to integrate 
our equations not along the real axis, but instead use a suitably chosen contour C such that the ingoing/outgoing 
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FIG. 1; The Priifer phase-functions (only real parts shown) resulting from a typical integration of eqns H24|l and H25^ with 
parameters I = m = 2,ujM — 0.5, a = 0.5M. Top panel: the integration of the P function was continued beyond the potential 
peak (located close to r, =0). Beyond this point its behaviour changes dramatically as a result of the Stokes phenomenon. As 
is clear from the graph, the derivative of P takes particularly large values at a series of points beyond r, — 0. This may signal 
the existence of poles of the phase function close to the real r« axis. That this behaviour could be expected, if one fails to 
account for the Stokes phenomenon, is clear from the equations given in Appendix 1 of ji^l . Bottom panel: as soon as the peak 
region is reached, we shift to the P function (the matching is performed at the point where the curve "breaks" in the figure) 
which is much better behaved, basically because it describes both ingoing and outgoing waves, see eqn. (1221 . Tests show that 
the result does not depend on the choice of matching point, as long as it remains close to the peak of the effective black-hole 
potential. 



wave solutions, eg. e^*"^""*, have comparable asymptotic magnitudes (of order unity) as r, — + oo. The simplest choice 
of integration contour is a straight hne with slope — Im lo/ Re to, see Fig|21 Similarly, a line with slope — Im k/ Re 
k will guarantee purely oscillatory behaviour near the event horizon. Such paths are par allel (for |r*| ~^ oo) to the 
so-called anti-Stokes lines which play a crucial role in WKB and phase-integral theory The anti-Stokes lines 

are curves along which / Q^/^dr* is purely real and as a consequence the WKB solution of ^ is purely oscillatory. 
It should be noted that a similar integration contour has been used in a calculation of relativistic stellar pulsation 
modes |49l |. 

We use a parameterization 

= r':'^ + pe'^ (32) 

with — oo < p < -foo and /3 given by, 

r._ ( -Arg{uj) for p>0 , .„„^ 
\ -Arg{k) for p < , ^^^> 

The matching point r™^ is chosen to lie on the real axis ( p = 0) and some experimentation shows that reliable results 
are obtained as long as the matching is performed near the peak of the curvature potential (the extremum of Q). 
Rewriting Eqns. (|24|l and H25|l in terms of p we have 

rIP 

^+e-*^(p2 + Q) =0 (34) 
dp 

^ + e-'^ (u; - ^] sin2(F + \uj\p + ujr^^) = (35) 
dp \ uj J 

The radial coordinate r will also be assumed complex, and we determine it from 

dp r'^ + a'^ ^ 
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FIG. 2: The integration path C in the complex r«-plane used to calculate QNMs. The slope /3 is chosen in such a way that 
divergent/decaying solutions at infinity and at the horizon are avoided. 



With the help of (O it is easy to see that 

r-«r+(l + p/'--) (37) 

where c = (r+ — r_)/2Mr+ as p ^ — oo, and that r « r, as p — *■ +oo. Hence, the proposed integration contour 
corresponds to a straight hne at large distances in the r-plane, and to a clockwise spiral as the horizon is approached. 
This strongly resembles the more complicated integration contour used in the phase-amplitude calculation of QNMs, 
see figures 1 and 2 in [l5l |. 

Not surprisingly, some care must be taken in implementing complex integration contours in black-hole problems. 
For example, we need to understand the analytic properties of the potential Q. From Eqn. Q we deduce that this 
quantity will have poles at r = ±ia. Since the function r»(r) is multivalued, there will be a series of poles in the 
left half of the r,-plane, symmetrically located with respect to the real axis. All of these poles have the same real 
part (which is zero for the non-rotating case), while their imaginary parts satisfy Im 3> Re r*. In practice, the 
presence of these singularities provides an upper limit to the complex-rotation angle. This is easiest illustrated by the 
Schwarzschild case. According to Q the point r = translates into 

n(0) = 2Aflog(-l) = ±2A/i7r(l + 2fc) with fc = 0,l,2,... (38) 

It is easy to understand that one should not rotate the integration contour beyond this angle. 

Another important issue that deserves some discussion concerns the boundary conditions in the complex plane. 
When we are solving the system of phase equations 13411 and H35|l we are actually imposing the QNM boundary 
condition on the contour C at p — > ±oo. However, the "true" boundary conditions are to be imposed on the real 
axis (r* —^ ±oo). It is perfectly natural to wonder if we are, in some way, altering the problem by using the contour 
C. This question has been discussed in the context of WKB and phase-integral studies (and the complex r-plane) 
|l8l l50| . It was concluded that it is safe to analytically continue the boundary condition in the proposed way. 



B. QNM excitation coefficients 

We have used the Priifcr transformation scheme described above to calculate both QNMs and Regge poles for scalar 
waves in the geometry of a rotating black hole. The method is easily extended to other classes of perturbing fields, the 
most interesting of which is likely to be that of gravitational perturbations. As will be discussed below the method 
is somewhat limited for extremel y ra pidly damped QNMs, and hence it cannot be used to investigate the asymptotic 
behaviour of the mode spectrum 27]. However, from an astrophysical point of view the slowly damped modes are the 
most important and for these our scheme works very well, indeed. 

In the case of QNMs one might also want to assess to what extent the modes are excited given a certain initial 
perturbation. As discussed in |40j the answer can be quantified in an initial-data independent way in terms of the 
"excitation coefficient" 

= ^^^^ (39) 

where we have assumed that A-^ ~ an{oj — ojn) in the vicinity of the n*^ mode frequency. While one can formulate 
expressions for this quantity within our previous scheme, we have found that the precision of the actual numerical 
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results is somewhat unsatisfaetory. To remedy this we have devised an ahcrnative scheme. The purpose of the 
fohowing discussion is to present this scheme, with the particular aim of achieving accurate results for the excitation 
coefficients. 

The first half of the calculation remains identical to the one discussed above, and we thus obtain a solution PL{r^'^) 
(say) which satisfies that boundary condition of ingoing waves at the event horizon. At spatial infinity we now obtain 
the appropriate solution in a similar way. Let us choose the phase of the square root in such a way that the solution 
corresponding to P+ « +iQ^/^(r*) has the desired outgoing-wave character. Once this solution is extended to the 
matching point we have a second solution P+(r™''). Thus a resonant overall solution, a QNM, follows from the simple 
condition 



(40) 



Now we want to determine also the excitation coefficients. We can do this by considering the solution for a general 
frequency, that contains both in- and outgoing wave at infinity, rather than a specific QNM solution. Assume (as 
before) that P+ « +i<5^/^(r*) corresponds to an outgoing wave solution, and that P_ w —iQ^/^{rt:) represents an 
ingoing wave. This means that, after integrating these two functions to the matching point the general solution 
corresponds to the linear combination 



(41) 











V'fl(r-r) = ^exp 


/ * P+dn 


+ Pexp 


/ * P-dn 



where the integration is initiated at r^. We can write the logarithmic derivative of this function as 



P+ + {B/A)P_ exp 


/;^^(P_-P+)d 




1 + (B/A) exp 


/;f (P_-P+)dr; 





(42) 



is that the logarithmic derivative is continuous, i.e. that the above 
expression is equal to PL(r™*). This condition can be written 



B 
A 



Pl 



■ exp 



Jj-OO 



(P+ - p_)dn 



for r* 



(43) 



Now recalling that P_ corresponds to the ingoing wave solution at infinity we know that a QNM must correspond to 
P = 0, and wc see that the above condition is the same as the one we derived previously. 

To determine the excitation coefficient requires a little more work. First we need to relate our amplitudes A and B 
to the asymptotic amplitudes Aout and Ain- Recall that the asymptotic amplitudes are determined from 



00 



(44) 



We want to compare this expression to our solutions. Since we have assumed that the WKB approximation is valid 
at the point where we initiate integration (rj°) we can use 



/ P±dr* « ±i / Q'^l^dr^ = ±ia;(r* -rT)±i 
beyond (the last integral should be easy to evaluate). We then have 



Q 



QV2 + ( 



dr. 



ip+ A exp 
And we can identify 

Similarly, we get 



/' 



P+dr^t 



= ^exp 



A exp 



j4i„ = B exp 



ia;(r* - r^) + i 



dr* 



for |r*| > \r° 



-iuir'^ + i 



(Q_ 



dr^ 



zwrr — I 



Q- 

9V2 



■ UJ 



dr^ 



(45) 



(46) 



(47) 



(48) 



10 



Collecting the above results, we have 
Ain P+ - Pl 



Pl - P- 



■ exp 



{P+ - P-)dn + 2iujr° 



2i 



(49) 



The final step in the derivation consists of taking the derivative of this expression with respect to to and evaluating it 
at a QNM frequency. This immediately leads to 



1 



Pl-P- 



duj 



exp 



{P+ - P-)dr^ + 2iur° 



2i 



(50) 



from which we can calculate the desired excitation coefficient. 

To conclude, let us summarize the additional steps that are required in a calculation of the excitation coefficient 
once a QNM frequency has been determined. First we need to evaluate the two integrals J^^ P±dr^. This is best 
done by introducing additional dependent variables, let us call them $±, and then integrating the equations 



(51) 



together with (|24|) from rj° to r™'^. Secondly, we need the derivatives and These quantities are also best 

calculated by introducing additional variables, i.e. by integrating 



dij du 



dr 



using appropriate initial data, i.e. 



at r^. Thirdly, we need to evaluate 



dio J 



dP± ^ i dQ 
'd7 2Q1/2 



dr^ 



(52) 



(53) 



(54) 



This is readily done by first expanding the integrand in inverse powers of r and then integrating. The result is 
convergent in all relevant cases. 



C. Numerical implementation 

Before concluding the description of the methods we have devised to study black hole resonances, it is worth 
commenting on the numerical implementation of the schemes. Overall we have a system of two coupled nonlinear 
first-order equations to solve numerically (for either choice of phase- function), namely, eqns H34() or (|35() together 
with 1)36(1 . In reality, since we are dealing with intrinsically complex quantities, we split each equation in its real and 
imaginary parts. At the end of the day, we arrive at four equations to be solved. This is done by means of a standard 
variable-step Addams method. 

To complete the solution of the problem, we combine the result of the integration of the phase functions with a com- 
plex root-finder routine (making use of Miiller's algorithm 51j) to search the complex frequency/ angular momentum 
plane for zeros of the ratio Ain/Aout- A prior knowledge of the approximate location of a resonance greatly speeds up 
the overall procedure. 



IV. NUMERICAL RESULTS 



The main aim of this investigation was to devise simple and efficient methods for computing black-hole resonances 
and scattering phase-shifts. As should be clear from the previous Sections, the phase-function schemes that we 
are discussing are straightforward to implement. They also turn out to be numerically robust and provide results 
of excellent accuracy. We have tested the methods by calculating known results for both Schwarzschild and Kerr 
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I 


n 




nM 


B„ 


2 





0.373672 


- 0.088962i 


0.126845 + 0.0206694 




1 


0.346711 


- 0.273915i 


0.048328 - 0. 2236184 




2 


0.301052 


- 0.478276i 


-0.190336 + 0.015118* 




3 


0.251504 


- 0.705148i 


0.080519 + 0.079964j 




4 


0.207514 


- 0.946845i 


-0.016727 - 0.060630* 




5 


0.169299 


- 1.195608i 


-0.001566 + 0.035783i 


3 





0.599443 


- 0.092703i 


-0.093672 - 0.049606* 




1 


0.582643 


- 0.281298i 


-0.152350 + 0.269067* 




2 


0.551683 


- 0.479093i 


0.414353 + 0.143038* 




3 


0.511956 


- 0.690338i 


-0.041146 - 0.412997* 




4 


0.470173 


- 0.915660i 


-0.220524 + 0.225274* 


4 





0.809178 


- 0.094164i 


0.064924 + 0.065661* 




1 


0.796631 


- 0.284334i 


0.263139 - 0.249790* 




2 


0.772710 


- 0.479904i 


-0.546231 - 0.439093* 




3 


0.739835 


- 0.683916i 


-0.323148 + 0.835511* 




4 


0.701524 


- 0.898240i 


0.867222 - 0.114054* 



TABLE I: Numerical results for the first few QNMs, and the associated excitation coefficients B„ , corresponding to gravitational 
perturbations of a Schwarzschild black hole. 

black holes. The obtained results serve as a validation of our numerical code(s). Having confirmed the reliability 
of our schemes we have extended the study to, in particular, the resonances of rotating black holes. Below, we will 
describe new results for both QNMs and Regge poles of Kerr black holes. These results shed interesting light on 
issues concerning the excitation of the long lived QNMs that are known to exist for near extreme Kerr black holes, 
the relation between the approximate QNMs for a ^ M used in (for example) , and also the effect of rotation on 
the black hole Regge poles. 

A. Quasinormal modes 

As a first application of our numerical methods, we have calculated the first QNMs of a Schwarzschild black hole. 
The obtained results agree perfectly with results in the literature , and confirm that our method is reliable as long 
as the modes are not damped too rapidly. The method breaks down when the imaginary part of the mode frequency 
is about one order of magnitude larger than the real part. This means that one should not expect to be able to use 
our scheme to study the asymptotic behaviour of the QNM spectrum. However, all astrophysically relevant QNMs 
are readily calculated. We have also tested the obtained excitation coefficients b y c omparing results obtained for 
gravitational perturbations (using a code based on the Sasaki-Nakamura equation j53|) to the phase-integral results 
in [s^ and the data obtained by Leaver A sample of results is given in Table By comparing the listed results 
to those in Table II of we conclude that the present results, which agree very well with those of Leaver apart from 
a few sign discrepancies already pointed out in |53 |. are reliable. In fact, we believe that the present method provides 
the (so far) most efficient tool for investigating the QNM excitation problem. 

As the black hole acquires rotation, each of the Schwarzschild QNMs (for a given I) splits into 2^ + 1 distinct 
modes. We have tested the reliability of our scheme by reproducing the Kerr QNM frequencies obtained by Leaver, 
cf. Tables 2 and 3 in T^. A sample of results for 1 = 2, obtained by implementing our Priifer phase function scheme 
for the Sasaki-Nakamura equation, are provided in Table IIV Al These results are in perfect agreement with those 
obtained by Leaver's continued fraction calculation for m = and m = ±1. Furthermore, Table is useful since 
the / = 2 QNMs may contribute significantly to the gravitational-wave signal from, for example, a black hole born 
following a supernova core collapse or a compact binary merger (see |m| for discussion) . 

Kerr QNMs are known to exhibit interesting behaviour as a varies |6j. In particular, we see from Table HVH that 
rotation has a pronounced effect on modes with m = I > 0. These can be thought of as "equatorial" and "co-rotating" 
with the black hole (cf. the symmetry properties of the associated angular functions, which in the case of scalar waves 
are spheroidal harmonics Sf^{9)). In the limit a ^ M the co-rotating modes have been found to become long-lived 
(Im uUnM — > 0) and accumulate at the frequency mu)+ ^, J4t.2.0j 21]. Hence, rotation has a significant effect on these 
modes. In contrast, the "counter-rotating" modes (m < 0) are much less affected by rotation and remain relatively 
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ajM 




m = 


= 2 




m = 




L 




m = 


= 




m 


= - 


1 




m 


= - 


2 










373672 - 





088962i 





373672 - 





0889621 





373672 - 





0889624 





373672 


-0 


0889624 





373672 


-0 


0889624 





10 





387018 - 





088706i 





380432 - 





088798i 





374032 - 





088898i 





367812 


-0 


0890044 





361768 


-0 


0891144 





20 





402145 - 





088311i 





388248 - 





088488i 





375124 - 





088700i 





362738 


-0 


0889354 





351053 


-0 


0891824 





30 





419527 - 





0877291 





397330 - 





0879951 





376985 - 





088353i 





358366 


-0 


0887634 





341333 


-0 


0891844 





40 





439842 - 





086882i 





407979 - 





0872571 





379682 - 





087827i 





354633 


-0 


0884844 





332458 


-0 


0891314 





50 





464123 - 





085639i 





420632 - 





086173i 





383318 - 





087069i 





351491 


-0 


0880914 





324307 


-0 


0890314 





60 





494045 - 





083766i 





435969 - 





0845641 





388054 - 





085995i 





348911 


-0 


0875664 





316784 


-0 


0888924 





70 





532599 - 





080794i 





455122 - 





0820851 





394129 - 





084453i 





346870 


-0 


0868844 





309808 


-0 


0887174 





80 





586015 - 





075626i 





480231 - 





0779551 





401917 - 





082156i 





345356 


-0 


0860034 





303313 


-0 


0885124 





90 





671634 - 





064863i 





516292 - 





069804i 





412004 - 





078483i 





344359 


-0 


0848654 





297244 


-0 


0882814 





99 





870861 - 





029510i 





572720 - 





0462091 





423685 - 





072701i 





343890 


-0 


0835504 





292107 


-0 


0880524 



TABLE II: QNM frequencies (lOnM) for gravitational perturbations of Kerr black holes with various spin rates. The data is 
for I = 2, and illustrates that the co-rotating modes m > become long lived as the black hole approaches extreme rotation. 
In contrast the counter-rotating modes m < are much less affected by the rotation. 

close to the initial Schwarzschild QNMs for all values of a. In Fig. O we illustrate this behaviour by showing the 
change in the first few scalar field QNM frequencies as the rotation rate of the black hole is varied. The data in the 
figure corresponds to < a/M < 0.99. 



-0.05 

-0.1 
-0.15 

-0.2 
-0.25 

-0.3 
-0.35 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

FIG. 3: The trajectory (Im wM as a function of Re ojM) of the first two I = 2 scalar field Kerr QNMs for varying black hole 
rotation < ajM < 0.99. The position of each mode at spin values a/M = 0, 0.1, 0.2. ..0.9, 0.99 is represented by dots. 

The behaviour of the QNMs as the black hole approaches extreme rotation is particularly intriguing. Unfortunately, 
our method breaks down beyond a « 0.999Af . The likely reason for this is the multivalucdness of the complex 
coordinate, cf. for example (|38|l . and the need to initiate the integration close to the horizon, see eqn H37I) . This 
becomes tricky as the event horizon and the inner Cauchy horizon come close together. A similar problem is known 
to arise in the phase- integral study of the QNMs of a Reissner-Nordstrom black hole [s^ . This means that we cannot 
push our calculation all the way to the extreme limit. Of course, we can determine QNMs for all rotation rates that 
are likely to be astrophysically relevant. Still, the behaviour in the extreme limit is interesting from a conceptual 
point of view. Especially since one can approximate the QNM spectrum for near extreme black holes |^. In figure 01 
we compare our numerical results for scalar field QNMs to approximate results for a « M and uj « mw+. The 
figure shows that the numerical results for the first two I — m — 2 QNMs, as extended from the Schwarzschild limit, 
approach the first two approximate modes as a —^ M. This is very encouraging since it demonstrates a correspondence 
between two rather different computational schemes. It also provides support for any argument concerning the nature 
and relevance of these long-lived modes that is based on the apgroximate formulae. The comparison also indicates 
the range of validity of Detweiler's approximate mode-condition 6] in a useful way. 

It has been argued that the existence of long-lived QNMs for rapidly rotating black holes could prove important 
for gravitational-wave astronomy. The basic idea is that a slowly damped mode will be easier to detect provided 
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-0.15 



1.01 



FIG. 4: Comparison of the numerical I — m — 2 long-lived modes (the first two are shown here, denoted by boxes) with 
modes derived by solving Detweiler's approximate mode-condition, (represented by lines, with the circles denoting the modes 
at a/M = 0.999,0.9995, ...,0.999999). It is clear that, for the fundamental mode, the two calculations are in agreement. The 
same seems to be true for the first overtone, although we have not been able to proceed numerically beyond a = 0.999Af in 
that case. The situation for the higher overtones also remain unclear. 



a/M 








m = 


+2 
















m = - 


2 
















uM 
























0.0 





483644 


-0 


096759i 





5252 


+ 





0950i 





483644 


-0 


096759i 





5252 


+ 





0950i 


0.1 





499482 


-0 


096666i 





6708 


+ 





2710i 





469295 


-0 


096708* 





4354 


+ 





0358* 


0.2 





517121 


-0 


096382i 





8355 


+ 





8432i 





456198 


-0 


096546* 





3790 


+ 





0132* 


0.3 





536979 


-0 


095839i 





8370 


+ 


1 


8608i 





444168 


-0 


096299* 





3412 







0044i 


0.4 





559647 


-0 


094931i 


-0 


9646 







0827i 





433059 


-0 


095984* 





3145 


-f 





0016i 


0.5 





585990 


-0 


093494i 


-0 


3913 







2772i 





422751 


-0 


095616* 





2950 


+ 





0018i 


0.6 





617364 


-0 


0912451 


-0 


1359 







2436i 





413146 


-0 


095204* 





2803 


+ 





0034i 


0.7 





656099 


-0 


087650i 





0005 







1699i 





404163 


-0 


094759* 





2691 


+ 





0057i 


0.8 





706823 


-0 


0815204 





0638 







0754i 





395734 


-0 


094285* 





2607 







0083i 


0.9 





781638 


-0 


069289i 





0432 


+ 





0143i 





387799 


-0 


093790* 





2550 


4- 





0108i 


0.99 





928026 


-0 


031069i 


-0 


0027 







0030i 





381041 


-0 


093330* 





2545 


+ 





0123* 



TABLE III: QNM frequencies (unM) and associated excitation coefficients (-Bn) for scalar perturbations of Kerr black holes 
with various spin rates. The data is for / = 2 and m = ±2. 



that it is initially excited to the same amplitude as a more rapidly damped mode. Since we now have the means 
to calculate the excitation coefficients for the Kerr QNMs we can analyse this problem quantitatively. As we have 
indicated in a previous paper [2l| , our results do not support the idea that the long-lived QNMs will be of particular 
importance. The reason for this can be seen from Table Ulll . where we provide results for the fundamental m = ±2 
QNMs for various values of the rotation parameter. From this data we see that while the absolute value of the 
excitation coefficient for the counter-rotating QNM only changes by about a factor of two as the black hole spins up 
towards the extreme limit, it drops by about two order of magnitude for the co-rotating mode. This indicates that 
the co-rotating long-lived QNMs may be significantly more difficult to excite in a dynamical process. 

We can use the results in Table Hill to discuss the "detectability" of the long-lived QNMs. Following we express 
the effective gravitational- wave amplitude associated with a single QNM (with complex frequency a;„) as 



, / Re ojn A°^' 
hcS (X X = 2 

Im CJn an 



Re ujr, 



Im u!, 



(55) 



Given our numerical data, we can compare the result for co- and counter-rotating (scalar field) modes. Figure 
(which is identical to Figure 1 in pll |) shows the outcome of this comparison. From the figure it is clear that, even 
though the co-rotating mode is much longer lived, its "detectability" decreases significantly as a — > M. We conclude 
that the present results suggest that the effective amplitude of the slowest damped QNM for a rapidly rotating black 
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hole may be several orders of magnitude smaller than the amplitude of the corresponding mode of a Schwarzschild 
black hole, or indeed the counter-rotating Kerr modes. The possible constructive interference of the many long-lived 
QNMs that exist for near extreme Kerr black holes, and which could still lead to a relevant gravitational- wave signal, 
has been discussed in [211. 




0.001 I ' ' ' ' ' ' ' ' ' 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
a/M 

FIG. 5: Assessment of the "detectability" of Kerr black hole QNMs as a — + M . The effective amplitude h^s (as estimated 
from Eq. is plotted as a function of the spin parameter a/M. We compare the slowest damped co-rotating QNM to the 

slowest damped counter-rotating one. The figure shows that the co-rotating mode (which becomes very long lived in the near 
extreme case) has a much smaller excitation coefficient than the counter-rotating mode as a — > M. 

As a final application of our method we have traced the path of a high overtone mode [n — Q,l — 2,m — 0) 
for varying black-hole spin. The result is shown in Fig. The purpose of this exercise was to investigate whether 
the QNM trajectory exhibits the characteristic "looping" behaviour, first found in the studies of gravitational and 
electromagnetic field QNMs for Kerr and Reissner-Nordstrom black holes [20, 56]. Fig. |B1 confirms that this behaviour 
is also present for scalar fields. It has been suggested [23,IE3| that this peculiar behaviour is somehow associated with 
the change in the potential near to the horizon as the spin varies. However, since it is unclear how the high overtone 
modes are related to the potential a fully satisfactory explanation is still missing (see, however, [s^l for an interesting 
point of view. 

-2.05 



-2.1 



-2.15 



-2.2 



-2.25 



0.14 0.16 0.18 0.2 0.22 0.24 0.26 0.28 0.3 

FIG. 6; The trajectory of the ninth I = 2,m = overtone for < a/M < 0.92 (we show Im uj„M as a function of Re ujnM). 
Dots represent the position of the mode for a/M = 0, 0.1, ...0.9. The characteristic "looping" motion is prominent for moderate 
to high spins. 
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In 


rn 


0.500 


2.086844 + 0.504735i 
2.188325 + 1. 4796131 
2.342582 + 2.392097i 


0.527566 + 0.019735i 
0.448169 - 0.888761i 
-0.350846 - 1.171626i 


1.000 


4.690075 + 0.501236i 
4.742988 + 1.494954i 
4.839323 + 2.466091i 


0.473472 + 0.575418i 
2.369138 - 1.190595i 
0.045487 - 5.147693i 


2.000 


9.889215 + 0.500325i 
9.915868 + 1.498744i 
9.967248 + 2.490034i 


-0.879437 + 0.579492i 
3.340145 + 6.574819i 
25.613332 - 5.812667j 



TABLE IV: The first three Schwarzschild Regge poles [In) and the associated residues (rn) for tliree different (real) frequencies 
ujM. 



a/M 




m 




+2 




m 






2 


0.0 


4 


690075 


+ 





5012354 


4 


690075 


+ 





501235i 


0.1 


4 


608103 


+ 





50071 H 


4 


763131 


+ 





500574i 


0.2 


4 


516055 


+ 





498888* 


4 


828131 


+ 





4987904 


0.3 


4 


412352 


+ 





4955824 


4 


885722 


+ 





495906i 


0.4 


4 


294773 


+ 





490504i 


4 


936387 


+ 





491909i 


0.5 


4 


160054 


+ 





4831944 


4 


980478 


+ 





4867484 


0.6 


4 


003103 


+ 





4728704 


5 


018235 


+ 





480332i 


0.7 


3 


815191 


+ 





4580884 


5 


049802 


+ 





47252H 


0.8 


3 


578931 


+ 





4356954 


5 


075234 


+ 





4631094 


0.9 


3 


247925 


+ 





3961774 


5 


094504 


+ 





4517944 


0.99 


2 


617843 


+ 





2781704 


5 


106494 


+ 





4396264 



TABLE V: Numerical results for the fundamental m — ±2 Kerr Regge poles {In) for a selection of black hole spins. The 
frequency is ujM = 1. 



B. Regge poles 

As mentioned in the Introduction, Regge pole calculations are comparatively easy to perform since as the associated 
wavefunction © remains well behaved as — > ±cxi. This means that one does not have to deal with divergences, and 
the identification of in- and outgoing wave solutions is straightforward. For a given complex Z, the set of equations 
(I24|l and H25I) can be directly integrated along the real r*-axis without any serious problems. 

To test our numerical code we have calculated the first few Regge poles for a Schwarzschild black hole for three 
different frequencies. The obtained results are listed in Table HVl As far as the Regge pole positions are concerned, 
the results are in agreement with the approximate phase- integral results given by Andersson [s^- However, our 
calculation has confirmed an earlier suggestion by Jensen "s^ that the phase-integral residues had a systematic error. 
As in the case of QNMs, we believe the present method is the most efficient and reliable way to calculate black-hole 
Regge poles. 

We next consider the Regge poles of a Kerr black hole. This is a problem that, as far as we are aware, has not been 
discussed before. As in the case of QNMs, the loss of spherical symmetry in the background spacetime unfolds the 
21 + 1 degeneracy of the non-rotating resonances. The pole trajectories (for the first two Regge poles) are illustrated 
in Fig. [3 for ujM = 1 and < a/M < 0.99, and a sample of numerical data is given in Table IVl Figure [7| shows 
obvious similarities with the QNMs trajectories in figure|3| This is not surprising given the close relationship between 
the two classes of scattering resonances [29^ . 

The most important conclusion we can draw from our results concerns the physical interpretation of the Kerr 
Regge poles. First of all, it is clear that m > (m < 0) Regge poles would describe co-rotating (counter-rotating) 
waves. Intuitively, we would expect that the "angular life" (the inverse of Im Z„, cf. "s^) of the former would be 
significantly increased due to rotational frame dragging. As in the case of QNMs, we find that the m — I poles are the 
most "sensitive" to the black hole's rotation. This strongly suggests that these poles can be associated with waves 
propagating around the black hole in the vicinity of the equatorial plane. More evidence to support this intuitive 
picture comes from the behaviour of Re In- The decrease (increase) of this quantity, as the black hole spins up, 
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is qualitatively expected from the corresponding decrease (increase) of the critical impact parameter and unstable 
photon orbit for prograde (retrograde) motion in Kerr spacetime. 



2 




FIG. 7: The first two Kerr Regge poles (we show Im l„ as a function of Re Z„) for various black hole rotation rates < a/M < 
0.99. Dots represent the position of each pole for a/M = 0, 0.2, 0.4, 0.99. The frequency is ujAi = 1. Note the distinct 
difference between the m > and m < branches which suggests that the former (latter) can be associated with prograde 
(retrograde) propagating waves. There are also obvious similarities in the behaviour of Regge poles and QNMs (compare this 
figure with Fig. I^l 



V. CONCLUDING REMARKS 



We have developed simple numerical methods for calculating black hole resonances (QNMs and Regge poles). The 
basic ingredients of our schemes are (i) the use of (Priifer) phase functions in place of the original wavefunction (ii) 
the integration of the resulting equations along rotated contours in the r* -plane (when the frequency is complex), 
and (iii) the implementation of Miiller's algorithm for locating the complex roots of 1/S{u!,l) function (the inverse 
scattering matrix element, whose zeros define the resonances). 

We have applied our schemes to the problem of calculating gravitational and massless scalar field perturbation 
QNMs and Regge poles for rotating and non-rotating black holes. The obtained numerical data verify all existing 
results in the literature with high precision. Moreover, we have provided some new results. In particular, we have 
shown how the slowest damped Kerr QNMs match the "long-lived" modes of a near extreme black hole which can 
be derived by an analytic approximation to the Teukolsky equation. This ties up some loose ends as far as the Kerr 
QNMs are concerned and also provides important support for using the approximate results to discuss the dynamics 
of near extreme black holes [21|. Our results for the excitation coefficients of Kerr QNMs add fuel to the debate 
concerning the detectability of gravitational waves associated with the long lived modes of a near extreme black hole. 
We have also provided the first results for Kerr Regge poles and discussed briefly how these results make sense given 
the established physical interpretation of their Schwarzschild counterparts. 

In conclusion, we find that our methods can be regarded as reliable and useful, apart from (i) for very rapidly 
damped QNMs, and (ii) in the case of near extreme rotations {a/M > 0.999 or so). This is, however, a small price to 
pay considering the profound simplicity and adaptability of the method. This is the main advantage of our approach. 
All other accurate schemes require a more sophisticated analysis: Leaver's continued fraction method jl^ is not easy 
to implement given the slow convergence of the relevant series, and the complex-coordinate phase-amplitude method 
developed by one of us [l5l | makes use of much more intricate integration contours than we used in the present analysis. 
Finally, we should note that, even though we considered mainly the case of scalar perturbations in this paper, our 
method can readily be used for the more realistic case of gravitational perturbations. We therefore believe that our 
"quick and dirty" approach to black hole resonances may prove useful in a variety of exciting contexts. 
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